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Ratiqmal Pade oonvetygwits to a Cimmtime =f Ube foreerd kW eure 
ator E are oamgited from coefficients @& its power series expences, 
These Pade oonvergents are used to generate weceruicn ferwalae aeiah 
are then employed to obtain aggefeaimete Humeriea] scletians 4f coventl 
second order differential equstions, A variety of integration goa 
lengths are used, Solution errc: is plotted against the Gulla af 
iteration steps and results are evaluated, It is «umillinied Ube Wcls 
procedure is beneficial as a method of deriving corrector faewmilas be 
be used in the predictor-corrector method of solving amd enter 
differential equations. A method for deriving ‘mproved predieter 
formulas is suggested. 
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1, Pivkerectuetaho. 

Eyal |3 has @aggweted a namerical aatied far mal ving aiffecw tie) 
eq@iation# «18 u@es rational Pade cawperpemts to appremiadle fumotLoee 
of the forward shift operator BE to generate reouesian fermilas. ii lw 
dmummistrated the merits of using such formilas in solving first oedar 
equatidma and has exhibited certain stability problems related therwte. 

The purpose of this paper is to extend the study of this metiuwl is 
the special category of second order differential equations in #aiam the 
first derivative does not appear, A computer program is deve]lim=l 
which generates a number of recursion formlas, applies then qumeesail¥e 
ly to solve the initial value problem y™+ y = 0, y() = 1, y’(0) = 0, 
and produces graphs of solution error vs number of iterations. The ba- 
sic program is then modified to solve the Coulomb Wave Byuation y+ 
4(1-2%- 51+) )20 with initial conditions appropriate to the selected 
values of the parameters 7 anc L. The Runge-Kutta method is used to 
compute eight values of the independent variable to start the solution. 
Finally, the results of these computations are discussed in deteil. 

Some of the advantages of this method over others in common use are} 

a. A rational convergent is generally a better approximation 
to a function than is the corresponding truncated series. This will be 
discusgzed further in the following section. 

be. Only about one half the total number of tabular differ- 
ences required by sethods using polynomial approximations are needed 
in this method for commarable accuracy. This is because differmques 
are distributed over both sides of the recursion forma. 

Cc. Since differencing is divided between both sides of the 











wrtiice, welleer alde ip (if Cereneel ee wey Cleese Tele comerey te 
meee of comtaliowe egCLrel Gt, HeyeLiT, Ge ope ay eevee, 

6. Seowecloy Cormalet CAP) Come Pete coerergeie che eee 
tommaly en@p bo @errecate Le co Giettal euupecier eevecrem, eke wil] be 
LLuwtrabet An Seotann &, 


2. Puirtivetition: of Bectratun Fareed. 

Duppeee it iw deelred to mlve Gre Apitial valee pepblen 
| yr(x) @ fie(x)), yx.) = Ye» y'(x,) = ae (#) 
antag Sins formally mi each wide 
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7 Pwr f Josce scan agar 
a he 


av get ei the left wide 
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Gabel’ pope - 4 fae) - ya) ° aor he ; 
Suisepilalig <, onde of variable 

(ex, + le,, di = hde, 

tec, dt = bee, 
are’ caf ining “he forward shift operator E gush that B f(x) = (aww), 
wm amt ww the zignt (tewating # ax an algvorsie grubiol)? 


Hh [ir she jine a) LD d2- Use aes & 4 idede 
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SettingEL]® 1 = 4. Se a eit te hip - eb Ce 


le 








AY La | 
\<- Ye : = (e-ih4, rhs, ‘ec . ») Je? ie Ve aie x te , 


fiona = mq? Oe 


In tetms of the backward difference operator v, def loal la vf) « 
iY, - Y;,, » this is 


9! how gf 
a “th me a ie =) CY i 


venir through by v ; 
h / — 4 s 
since W operating on a constant yields sero. 





For simplicity of ow we let a vy; and write 


2.2 $: 


teat a 
_- 
Let a.+ av + av + ose be the power sextes expansian of fl 
The power series can be approximated by a rational frastign hp wee ef 


the following technique. a + av + av + coe HAV er) etree 4 
leeVeray"s + cp® 
Multiplying each side by the denominator of the above freetiai an’ 


equating coefficients of the like power terms in v we pemeruie ptq+i 
linear equations in the ptg+l unknowns D, eee a> CO ceeSoe 

b= a, 

b, = a, + a, b 


°C 
b,= a,+ a,c + a, 
b= e,+ ac +a 


c, 
c.+ ac 


a eae > 4s 
The solution of this grstem gives us the ooafficionts of the retin 
al fraction which agreee with the series ost to P+q+!| term. 


It mould be noted that a given polpummlal of Cegree r cam be equate 








tieitioslily to ee@i of r+ 1 dimiiant fraviioee. Trewtetline Vie cals 
nowtiel afler 1, 2, ..., n vari, Wuowsivel>, @e cPdeda « collewtiae 
of polynomial fradtions whisk San Se arranged as « Val Mais of vuttem 
al convergents. Fig. 1 is the Pade table asscniahe! wis ey : 
where a third degree polynomial is used to paemate te Petinnal 
convergar.ts. 

The power series expansion of ; wi Fl , ree Ge sevet 
terms, 1s the sixth degree polynemial 

jsve tyra ivs tilg* +7 Vv o + Bog” et) 

The coefficients of this polynomial were used as imgete vs Pregeee 
& UCSD PADE (Naval Postgraduate Sohool Library Routine) to a@rbeis 
coefficients of the 26 Pade convergents whose numeratuts ex! danawifeters 
are polynomials of degree q and p, respectively, and mhere prq ¢ 6. 
Fig. 2 lists these coefficients, rounded to three signifiemnt figures. 

As stated above, each of these oonvergents is identically eqaal to 
(3) out to ptqtl terms, Moreover, a given convergent expanded algeutwaie. 
ally into a power series generally continues to agree qrite well with 
the related infinite series for a number of terns beyamd the (prqti)™ 
1.@., a rational convergent generally is a better approximation to a 
function than is the truncated series from which the cayvergent wus 
derived, A 

Dosigpatp a Sa mal rer ‘gal te iegcicr? 


AN! )y (te) 
ae rye met ys — 


Suppose we wish to find an ayprogiimete numerical eeletian of the equation 


by 





os ee 
y" «@ f(x,y). Sstitoting (4) for mt Uae in equation (2) we bare 
s « n't /| bul ts bales! 


Spcat + Ae (pe? 


” -. 4 P 
m——Fitties, sie (ecceesg iA Ln Ce i | Sarin amwult’ Mi 











hike Kel dh = OS, (oe nen eepreee Lhe Gre eye 
savrpitly, +t) Ea» WC, fe f *) nef vit he eableath 2 (4) 
Comte Che woul Taeiewhe Any (5) a@ Follow: 
ats) > () jae are) 
ra) BC G41) 
fubetitaling ‘er cowfficimets in (5) we hee: 
Q(Gje - Q(iley + ..-f aa, b'[ PCO), - a DY POP) =, 
7 — atte ~ ooo HR + n*[P(o)s, Ste See Ji (8) 
here e, be given in Lemme of coefficimts of tre Pete cope argent 
and of previously ommputed values of 2 and f, Recalling that 2,= vy, * 
Y,~ 2,,+ Y.. ve my now calculate y,#2,+ 27 - y.,. The value of 
f, uwwt in (6) will have bean onmm—uted from a first ayprve@esticsa to 
¥, emverated by a predictor formule. 


3. Stability Considerations. 


The recursion formula generated froe a given Padé comverfent is af 
the foray Q(0)a, = Q(1)a = of Made, + n'[P(O)f = ... tP(pIt. , } ee 
wnene Q(i) and P(4) are as defined in Section 2. 

Aubetituling vy, * y,- 29,,+ y), forain (7) we obtain 

Q(0) (yy BR, t¥0 9 ) # AL) Ce, 2H, Ag) ~ vee FUG) are 
Pig, *oy) * h’{P(O)f,~ P(L)fint vee? P(p)f.-+| 
Multiclying ere] collecting tems, this benamas 
alaly, * Farr) + co] x. = fece) + mata) + Q(0)] rater 

+ (<4) | Cg) + M1) + a(q-2)] n+ (=t)F [2d¢~) + et) a 
+ (at) ata dyy [CG ~ --e4 (-1)" PC )E1, | (o) 
The yemeral adetion of thts 4ttrormme eqieticn ip of We fare 


5 





awe Rese ~-« * of, ere Fon... Be 
are rowte of the ewem lated dareéterietic azreniic: ((e b Intrniered 
hater) emi the care counmtants. 

Conpider the 2ad order Pade caWergeet ne a (Forme 1) 
in Fig, |). From the digmussion of Section 2 we kiow tet tie tumeiae 
formula gmnerated by this campergumt will have the follewing soul iajaetey 

Q(O) = 13 Q(1) = 0, 121 
p(0) =43 P(1) =~ P(2) #2; P(4) = 0, 1 #3 

The general recursion formula (8) then becomes 
Y= WHi- Kat lat +75, titra} @) 

In the initial value problem (1) it ie seen that f,= «y,. 
Substituting in (9) and collecting terms: 

(1 +ib dy; = (2-28 dy (14Gb Dy, 
or, 
(12¢h )y.— (2-10 )y,,+ (12¢h )y,,= 0 


Finding roots: rae a4 -soh2\ t V Cae -wh*) -4(129h") 


aC rmrar*) 


- (1a-5h?) + ay \4¢4-241* 


Deh 


(10) 


Ignoring temas of the order of h’ we have 
vie tee - ,+2 

a = 
The general solution of (9) is of the form 


(11) 


y.Yorn+o,r, 
Ww determine the ocefficients oc; from Uke soletion y = oes x af the 


initial value problea., 


4 *#ert+ore@l (initial comticiews give y,* 1) 
y @or+ ora tat (y emma b © ised)... . Few aml) 
2 Wuizes af b,-vSe a goad onptet ae len) 








Tees wo have; 
ge = tk fusbty ge Zee TEAST, (42) 
Tils can be writt# in trigescmiric form 
se PELE) (Creer eS] GENT) oe ew] 
h 


where SS: One Bw 


~~ ea . 
yn: (Vier ) (4 + Ba)(cos woe rsinme) + (4-$1) (cos mO-1 Sun no)f 





= (Yr) (cos wo-h sin 78 ) (13) 
for mail h, hr® in wialdh pase 
an? (Vary (com hese. oh ) 
and this solution carwerges to the exact solution, 

It can bo wean from the graphs that Formulas 10, 11, 12, 15, 2%, 
oul 25 are the moot ptestising fove the standpoint of miniems error 
amplification. Unfortemstely, their sessociated characteristic equstians 
(fom. “rou eeventh degree) are too intrastable to permit analysis, 
It wuld be quite alaple, Jowever, to rbtein the roots of these 
Pyusetiong with an appropriety wompater program. The erm siaracteretics 
of Ue ayustione could Shen be determined fromm a siuely of the prepertion 
of tveir reece. Tiis is diatusmmel in Bettion 7 ap a project reoummembaf 
fer fdture etary, 

Tae foreinleg pa? eg teats Ow coeeerped wild Lhe eccueeletion of 
Lraeeion acoer withou) regerd for the fect tha eepece Le y,, an 7, 2 
fe (3) #1 affect cibeenees! armies. Te sancifer thle effort a 
differen \ ameceach «lll be tees, 





in beret of the U(Tereene eyushien (7), 
i] = . a — ‘ 
ies Cobact CELN ie eer en) a (ibe 


worse rte Une value of pclwtten of Gwe Uniti) welce creslee of ©, | 


oa FT. = Cle, Fi, TCs) 40 She Coveetion eftur inirates ce dellee 


, ff ‘| ‘| | - : 
_— a approciests | —— \ 2/5) Le tee powetleEf sevar. 


We CaPine Efr,1) © Uh) + O41). 

Cwfiniew 6, % y. - y,, We mabtract (9) frm (14) aed cytadn 
Gis PLeemdes +] SC ay EE) ot OF 2. | sb oe 

Aweuming tHet £)(x,,y)is differentiable with respect to y in 
the intmrval (¥,,y,) wi mg use the @man value theorms to write; 

(iect:)- 4, (Ks 8) (5i- 44) = én f, (15,8,) , 
if f (x7, ,) does not vary greatly we may let f, (x>,1;) = K and write 
€é,: 16, -€, m th?kK (ob €. + 12 €,,, ee Ca. | + & Choa) : 

Collecting terms: (i-+:h*h)é. = (as tee) ey. - (1 LAK) ey + EG) 
OF, (Meh ke, (2H se heey, = (1a-h*e IE... FECA). 

In the solution of the initial value problem (1), K = f,(x,,ti) = <2. 
Hence, (12+h’ )¢,= (24-10h' )e,, = (12th’ )e,,+ E(h,4). (15) 

Under the aswumption that E(h,i) is relatively insensitive with 
respect tc i, we may proceed as we did in solving (9). (As applied to 
this oase the trimcation error T(h) = ch’ y‘‘te,) » and y (#,) varies 
slowly when the h increments are small), We will hen@eforth asmanm 
E(h,1) to be a omastant: Dh’, 

Wemmst solve (12 +h )éi = 20h’ )é:.- (12¢h’ )é.4 Dh® 

6. = 0 (no error initially), <ul 
ec (assumed to be non-sero) 
The solution camsists of the guneral solution corregMaiim to the 


humigeneous part ¢ = co, r+ o,r” ml a particular smiisticn Gliese. 











porcdeng ao Une Ag- emrgeeues part s. = perieciar @letion ee be 
toecd Sy seweelre Great’, ie o constant &, = hel 
(ideh? A= (M108 )a - (lee )A + Oe. ee, 


a | 
A - a —— Dh = les 2 64° r 
= ( ob") af OG -m ) ra” r2 
The Complete selution i# thangs or+an +5 h. (v0) 


Let r,™ a + ib, r= a - 4b, 


Dendte c =i + if, Clad - ip, where d= —_ Sana As beset | : 


The solution (16) is transformed into trigdénasetric form tr letting 
Q@ are sin __ 


Then, € = |ri (24cos n@ - 2fain n6) “2 h*. Expanding: 
& o , 
é,=/rl b-2 bh cos n@ +i (S+ (a-1)2 h*) sinn® +Sh. 


In particular, if we spproximate 
a=i bh, i.e. r + i+ dh, r, = 1 - th : 

<. = (+h re -2 h cos n@ +—(f adn no) § “i h- (17) 

where @ = aro sin 7tr - For small h 
6. ¥ (fan y }-2 noon nh + Lain nh § ne h (18) 
This solution of the error propagation equation qualitatively ex- 

plains the behavior of the experimmutal remilts displayed in Fig. 16. 
a. If h is small, the sine term of | = 2 ni o0s nh + Len wn} 
2@ @ignificant and the coaine tarp is negligible. Etrur would be 
expected to vary nearly sinuwiidally. | 

b. 2, sboald be avrmcimetely equal to cerw adhere ah i¢ eal- 
tiple af¥, el:co the teme/ ain Gh ill venieb ot Khmer pelts, dines 
he 02 ant eleing for n © ra » kw, 2, ..., S, wo mmpmet be fied 








reso—eect etter woerw 6 = 147) 544, Ol, 62 ced TES. Pig, 1 we Cede 
bc he the come. Dror celle aes to Figs. ¥) ae bi om fe eeplained 
[imileris be eleier for n #105 h @ .0) Gh .05, rompemEively. 

c. ‘Thw branketed temp fy (14) is saltigliet Griel> wees, 
the ermor eacpitede emeclé grow with insiwemiae n iff) 1 eet hele 
ifpricl, The melificathas factor (7TH) in (10) is aly en append. 
@ation of iri. The eetual roots were femme to be 


at Ge p | Wiee- sen? 


saah* 


Inj = (ae-iteh™+ ash a h*(jan-Soh") 2 4 


1¢*#+ Ie Kr, bh” 

for all h,. Wepever h is small enough to insure that (16) closely agigumie 
imates (17), error maxima should remain nearly constant. Fig. 33 mie 
etantiates this fer h = .01. Increasing error maxima in Figs. 19 amd 41 
(h = .02 and .05, respectively) are attributed to the failure of h to meet 
the criterion statmi above, 

Yo explanation is offered for the high frequency error oscillatia} 
so prevalemt in eam of the graphs, e.g. Fig. 16. Investigation @ this 
pheqametor is muceeeted as a project for future study. 


uh, Demriptae of the * Pro 





This progr t¢ deeigned to solve a given second order differen- 
tial equation %h tits, suscesively using each of 26 Padé canrergmits 
to pwnerat* a central regursion formula. Integration of the ¢ived 
equellom ie te the preilctor-corrector method, using the Milae preeilaler 
formmla; oa Lala Ge Jae oh [557 ya | fecessary 
Wiarkire values acy caormcoted bp Une fucyee-Kutioa method. Fhe iterative 
prccwen ear he cacred oat so aw emer aw 900 ebece, atlaiton errar kaing 


40 





Seheeueet at cali sawp if am aepet gelethce of the powen epeecken be 
kneen. (raphe are plotted with error as setiieete aul came =f i beretdene 
as abscie®#a. Fig. 3 explains tie GGBéls yams in Ge pregreee. Fig, © 
is a printout of the progres cfd to evalaete re aqeeiies gereretlng 

the cosine fmeation, y* + y = 0 (firegial TNO 2). Fie. 5 is «a price 
of the program used to evaluate the Coulum Wave Gyuthen 51, (/-“ joa 
(Progres TMOD 6). 


5. Dipgwamon of Data. 


a. Program TMOD 2 

(1) Figs. 7 through 25 are graphs of progressive arr Pegilt- 
ing when the initial value probles y" + y = 0, y(0) #1, y*(0) = 0 is 
solved using an integration step length of .02 and corrector fevmiles 
derived from the various Padé convergents. Fig. 6 relates the recureian 
formula number shown on each graph to a description of the correqpemdiing 
Pade fraction. Formulas 16, 17, 20, 21, 23, 24 and 26 diverged rumiidly 
and no associated graphs were plotted. 

(a) Formulas 10, 11, 12, 15, 22 and 25 preduced remarkable 
accuracy over the entire range of 900 iterations. 

(b) Formula 6 (Fig. 12), although less acwwrate Laas the 
group just distussed, has an iapertant property: error maxima are 
d4minished with increasing iteration. 

(c) The remaining graphs indicate that error erie are 
not being dimeed bat are insreasing linearly with reapeet ton, the 
nemiper of iterstione, 

(2) The creek ah of data was obtained by repeating the pow 
cmtioe sccqoerieuey: af ier deorpeseing “ve interreticon step lemgt) to .91. 
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The following wew sbearvel) 

(a) Taw cane oiler proved to be cot somes BLE « Sep 
deewty of .0c (Mwreked AO, 11, (1, 14, Sh aot 75) one Com oo ne 
Lees eootrate With the emailer wher leneth (centrary Go sil eould be 
aipecte’). Purnilas if, £2 a9 2§ shew a onal] ofewr eeyplifieation while 
all the othery in Chie «ree Wow a. cewh Coaster errer amplification. 
Fige, © tare 33 «re eraphe owing Uiems remlte. 

(b) Pereale 6 provides le@e acctirate regults thas with the 
@restwe etep leneth bet ite error damping charecterivtics mummies. (Fag. 27) 

(c) Forwila 5 is are sauble. Its asguracy is compirible 
To that in the preceding own. (Fig. 26) 

(d) All other formulas proved to be @ore accurate Liay sti 
used with the greeter st#p lesg@th. All of than, howewer, amplify orrers. 

(e) The @rrer granis a@goeiated with Formiilas 9, 10, 11, 
13, 14, 18, 22 amd 75 pere idesitical. (Fig. 33) 

(3) Increasing the step lergth to .05 had the following effeots: 

(2. Th® formulas whisk proved best with a step leagsh of 
.02 are mili lees ecctirete with @ step lumgth of .05. Figs. 34 Sliecugh, 
1 are the error grwohs asemcleted witS Formilaz 6, 10, 11, 12, 13, 15, 
22 and 25 weieg the .05 wlep lemgth. These grajis Miow the only forutilas 
cwlaining eey wreve daeolnc cllaracterigtics to be nigiters 11, 12 Gm 25. 

b. Pregraa Damp 6 
(1) Tie Comlmab were Gyustion we» eilved for two compinaticns 

i wolese of Leet 7 (lel, 61; a Let, 7)), eet a Weriond Lote 
eraticn slep lesetoe. IX wen col feeeiale tv eenerete error creche far 
Leeot eyustioce. Fle, Of Aece eaccinete ef devlabion of competed sel. 
toon Looe (owe wuss aber We nomter of iterations tnfieated, Das 


4? 








“) Uimdtatieed ar ths tables wet} 3], the Winimee esses Alececpetle wad 
1A (le Tom ). A> @@ewivetiaw of Fir, G! wevealo She folboueg: 

(a) Identical re@ilts Mare citalned 4) all Gates ere 
Furmies 9, 13, 14, 15 amd 18 were vem. Dis grou elem prumiuesd ‘ike 
wuet accurate solutions. Jlieduming step length Lmprove’ eoqureay Ukewagh- 
out the range. 


(>) Formula 8 solutions wee slightly less aécumate ‘than 
thowe of the group above in solving the equation y" + y (1-3-4) =z 0. 
Solutions of the equation y" + y(1-— -<) = 0, using a step length 
of .02 were more accurate than with any other formmls. 


6. Conclugjora. 


In his Ntmerical Methods[3|, Kopal illustrated the shortcomings of 
going-ahead recursion formas; those derived from the relationship 
i ad fr. - When uaing GolmeregEpee formles stemming from the 
integration formula y, = +f f;, an spproximation of f, must be ob- 


og@ 
tained before y, can be Aiea Therefore, the method suggested by 
Kopal is essentially the predictor-corrector method coupled with a pro- 
cedure for deriving an wmlimited supply of corrector formulas. Some of 
these we have seen produce quite accurate solutions to several differen-~ 
tial equ@tions of second order. A glance at Figs. 43 ~ 45, the error 
of y" + y © 0 after solution by the Milne method, is evidence of the ayaa 
jority =f a neoier of the correctors derived from Pade convergent. 


It is mumgested at this point that the less accurate geing-ahead 


tyre foateulte aiyot be g@merated in the Mame way and Gsed as predictors. 
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Toe fait Lt Sertain sonrectete eriter Dees Made aererreyis ert ton 
ety of pretieter 6erivetian., 

Vom f Ue correotce fcommled dieweses! ip (244 pepe ase wubdle 
for 215 «jettlons over CL Goa oP oh. befitiees! fomeedee denid be 
darko! from « sclmemal <oprontaniitoy ta) 2)" op tse tegee, 
tw Take ama dees not appear promheliag . In dealing will asemel ated 
OharsaLeristic eywakiove of hieeer depree it will be ectremily UAs¥iaul4 
to locate valuew of pase bers e@he8) will retails all roots wileia the 

The meuwliest reslts obtained with Pewmele 15 (Fig. 21) Sadianke 
thet o frumth Gecree polnestal areremnaation to fC sy trian 
to iownire “het truncaticn eFror i mut a @igefieant peoblen, Th eral 
coe, Cherefore, tiat Che ome of polynomials of higher decree ewald ims 
vite acditiunal e#tebllite camplisatians. 

In wommury, it ie conmmluded that the procedure suggested ty Eapal 
shoald be viewed of a method of deriving predictor and osrmeter farmalss, 
thee remilce eashievred with see of Che correctors derived ts tis Sauer 
Aodioate thet thie prosedupe ie primtiaing Ant merits further stualy. 


?. Bieommeecia tions fot Pybire siudles, 

». LDvvretiqnte the roole of the darecterieice sqeetians asmeiaial 
wish whe mew proeleing cecucebee farmilag Glecussed in this paper. Ae 
vertetn eel cabclate veluem of & for wnian all voete lie within the aut 
virele. Wit) tet dete evaileule, ene soci neatly selsat the bert 
reeetrusn, frerecis for solving * poe differeitias equgiicn ower a 











pom tel warnee of ealuee of Ue beclepereten! vari ele, : 

». Omirlre « Pade table of euliendl eererercert« =} tit ea te® | 
free which te coteln « nec af ewinpediee!l reerreion (oemion, See 
wt vuhe each of Cores vo tae fer tiie prefieter fort. cual is Fee 
PO >. Calcel@te errors ay) cog with Pie. 42. Dh thle some, 
evaluste the merite of predictor frimmilet derived froe ration) ee 
cOnverceits. Ralelon 6} Cemeribe# tire adi ar tare of ming preadistere 
eno correctoré of the seme order, 

c. Extmnd the polynomial sepremmitation of {" herent 
ewren terme aid derive higher order recursion aemse Inrewtigate 
the roots of these formulas as miggested jn subparegrégt a aleve. 


j= aul ie 
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wollte Pevee of (xleee Wf Ce Jiteeencent septate. 

b. Decker & Pact Gale of cwUleesl sovteepente taf (1e)”* 
Creat wiles to shtear © number @C peieeedieeld cocerstae Peeeiies, Bee 
@titute each of thane in turw fer the peeiictce Cormul, ued iw Pregeee 
TaD 6. CAlégdLite errors ami compare wit Pig. 52. Ty th9 eames, 
wveloets te sertts of jemdictor formulas Gerivet fhe notional bade 
canvelrg@aets. Raleten | dewerihes the eivuntacw of axing aredictote 
and correctors of the same order. P 

co. Extemd the polynaeal represeptation of : — ; ene! 
@even terme and derive higher order recursion formilas. Lvwestigute 
the roots of these formvlas as @iggested in subparagrag « é50Ve. 


i. 
a= i 


™5 





Tir. 1 - Geepic Padé Table 








The {wt “het Urres of u 2e omnes are 
jul ip Snioeifjeste), 











4+ 


Vc, 


1 


- 


too 


co DB NN DH we 


10 
11 
12 
13 
14 
15 
16 
L7 
18 
19 
20 
21 
22 


23 


25 


. heed ® 
- Co@T4elervte of Rese Camverpenta Le ; ts : 


fase tor 
1,0 
1,0,1.0 
1,0,1.0,.917 
1,0,1.0, .917, .833 
1,0,1.0,.917, .833, .761 
1.0,1.0, .917, .833, .761,.700 
1,0 
1,0,.083 
1.0,.091, ,008 
1,0, .087, .003,~. 004 
1,0,.080,~,003,~-.010,-.005 
1,0,,074,~.009,~.016,~.011,-.005 
1.0 
1.0,0.0 
1,0,.600, .050 
1,0,-1.29,-,121,~.014 
1,0, ~.830,~.082, ~. 013, -.002 
1.0 
1.0, .083 
1,0,-1.0,.050 
1,0,-.954, .077,.002 
1.0 
1.0,<1.0 
1,0,-1.0,.123 
£50 
1,0, -. 877 
17 


> 
vt Let 
1.0 


1.0 
gO, <220 

1.0,-.917 

1,0,=.909 

1.0, =.913 

1,0, =.920 

1,0,-.926 

1,0,-1.0,.083 

1,0,-1.0,.083 

1,0,=.400,-. 500 

1,0, =2.29,1,25 

1,0,-1.83, ,832 

1,0,-1,0,.083 
1,0,+.917, .063, . 004 
1,0,-2.0,1.13,-,133 
1,0,01.95,1,12,-.155 
1,0,-1.0,,083,0.0,-. 004 
1.0,=2,0,1.068, -. 093, -.00% 
1,0,-8.0,1.21,-.908, em 
1,0,-1.0, .083,0.0, - 004, ~, Oo 
1.0,-1.08, .960,-.073,-, 00%, -. O01 





fy. 3 - Legené of Program PPtele 


at | 
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M 


IN(I) 


ID(T) 
DEL(K,.,J) 


QO0 
W(J) 
ro 


P(d) 


lnhwomyient variable 
Dep~myaent variable 
y p 


y" 

y 
y-y 
vy 


wy 
Coefficient of 2 numerator term in Pade onnrweary amt 
Coefficient of i denominator term in Pade (ory orwpen t 


Initial values 

Number of first convergent to be evaluated (isupal) 

Number of last convergent to be evaluated (162g) 

h (AIntegration step length) 

Mweber of starting values to te computed by Mimge-Liitte site? 
Total number of iterations desired 


Denctes number of terms in nuemurator of formia mo. (I). 
Used for information in printout 


Seme as IN(I), except refers to denceimster torus 
Array which stores Pade coerficients. 
nemerator coefficient; Kw2, coef *iolemt. 

I is intlex no. of convergent. J is the numer of a 
mecific tar. E.g., 0HL(2,7,4) dumgtes the caf fialeit 
of the fourth tems in the denemineler af comperpont ae.7. 
Comfficient af #, in recumaien formule 

in rewurakle: formals 

ie mwrureden formle 


G=1 denstes 


Com@ficient of «,.- 
Caaf icient of 9 
CuiPTieliest of yr. fa munwncies favemla 


ieee os. SS sovrerpeet ew its aeecctesed receretian femme 
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1 SeSelesLlAST) 

If (L=-IQ) 24525525 
24 Laet+]) 

GE TO 26 
25 PRINT 99% 


END 
FUNCTION MIMOM (MAN) 
AMEM 


IF (M=0)9015.3%25302 
DO? IF (N-1)498,804, 10% 
305 IF (MN) 3'%1 903536 
#96 CON? JMUE 

K=h-1 
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